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ABSTRACT 

We observed mid-infrared (7.5-22 /im) spectra of AGB stars in the globular 
cluster 47 Tuc with the Spitzer telescope and find significant dust features of vari- 
ous types. Comparison of the characteristics of the dust spectra with the location 
of the stars in a logP-K-diagram shows that dust mineralogy and position on the 
AGB are related. A 13 /xm feature is seen in spectra of low luminosity AGB stars. 
More luminous AGB stars show a broad feature at 11.5 /im. The spectra of the 
most luminous stars are dominated by the amorphous silicate bending vibration 
centered at 9.7 /xm. For 47 Tuc AGB stars, we conclude that early on the AGB 
dust consisting primarily of Mg-, Al- and Fe oxides is formed. With further AGB 
evolution amorphous silicates become the dominant species. 

Subject headings: dust — Spitzer — stars: AGB and post-AGB — stars: mass 
loss 



1. Introduction 



Stars on the Asymptotic Giant Branch (AGB) represent the final evolutionary stage 
of low- to intermediate mass stars. The large AGB mass-loss rates play a key role in the 
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cosmic circuit of matter. This mass loss consists of both gas and (sub-)micron-sized sohd 
"dust" particles. Dust in the circumstellar shells of AGB stars leads to a significant change 
in the overall spectral energy distribution compared to dust-free objects. It attenuates stellar 
radiation in the blue and visual range and re-radiates corresponding emission at mid-infrared 
wavelengths ('infrared excess'). The spectrum of the infrared excess contains broad spectral 
features characteristic of specific dust species. 

Up to now, dust properties of AGB stars could be determined only for samples in the 
galactic disk and bulge (e.g. Omont et al. 2003). The IRAS and ISO missions created large 
databases of mid- and far-IR spectra of dust-producing stars. Analysis of these spectra led to 
basic insights e.g. into the relation between dust properties and stellar C vs. O abundances 
or total mass loss rates (e.g. Molster & Waters 2003). However, as the disk and bulge 
samples are inhomogeneous with respect to individual metallicities and masses - which both 
influence the production of dust - several questions remained open. One main open issue 
concerns the interrelation between the dust composition, stellar variability, and evolutionary 
status of an AGB star. In contrast to stellar samples in the galactic disk and bulge, globular 
clusters (GCs) make possible observation of stellar samples with well-defined, homogeneous 
parameters. 



2. Sample selection and data reduction 

Among the GCs, 47 Tuc (NGC 104, [Fe/H]= -0.76, Briley et al. 1995) is especially 
suited for a study correlating AGB stellar parameters and dust spectral signatures. This 
cluster has a well populated AGB with indications of dust formation due to considerable 
IR excesses (van Loon et al. 2006, Ramdani & Jorissen 2001). Lebzelter et al. (2005) and 
Lebzelter & Wood (2005) were able to assign fundamental or overtone pulsation modes to 
most of the 47 Tuc AGB variables. Overtone pulsators exclusively populate the lower part 
of the AGB, while fundamental mode pulsators are only found at high luminosities. The 
NIR velocity amplitude was found to strongly increase during the transition from overtone 
to fundamental mode and to further increase along the fundamental mode sequence with 
luminosity. Based on these findings, we selected 10 targets distributed over the AGB pul- 
sation modes and luminosity range for observation with Spitzer. One star, Vljjl], shows a 
short primary and a long secondary period (Lebzelter et al. 2005). To this sample we added 
V18, a star located between the fundamental and the overtone mode pulsators with a strong 
infrared excess (Ramdani & Jorissen 2001). It was suggested by Lebzelter et al. (2005) that 



Throughout this paper, variables are named according to Clement et al (2001). 
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this star is currently in the luminosity minimum following a thermal pulse. The upper left 
panel of Fig. 1 gives an overview on the location of the various variables in the \ogP-K 
diagram. 

For these 11 targets low resolution Spitzer spectra have been obtained between 7.6 — 21.7 
fim during Cycle 1. Our spectra are based on the droopres products processed through 
the S13. 2.0 version of the Spitzer data pipeline. Partially based on the SMART software 



package (IHigdon et al.ll2004l ). these intermediate data products were further processed using 
spectral extraction tools developed for the "Formation and Evolution of Planetary Systems" 
(FEPS) Spitzer science legacy team. The spectra were extracted using a 6.0 pixel and 
5.0 pixel fixed-width aperture in the spatial dimension for the observations with the first 
order of the short- (7.6 — 14 /im) and the second and third order of the long-wavelength 
(14 — 21.7 /im) modules, respectively. The background was subtracted using associated pairs 
of imaged spectra from the two nodded positions along the slit, also eliminating stray light 
contamination and anomalous dark currents. Pixels flagged by the data pipeline as being 
"bad" were replaced with a value interpolated from an 8 pixel perimeter surrounding the 
errant pixel. The low-level fringing in the low-res olution spectra at wavele ngths longer than 



20 //m was removed using the irsf inge package (ILahuis fc Boogertll2003l ). The spectra are 
calibrated using a spectral response function derived from IRS spectra and Cohen stellar 
models for a suite of calibrators provided by the Spitzer Science Centre. To remove any 
effect of pointing offsets, we matched orders based on the wavelength dependend point spread 
function of the IRS instrument, correcting for possible flux losses. The relative errors between 
spectral points within one order are dominated by the noise on each individual point and 
not by the calibration. We estimate a relative flux calibration across an order of ~ 2 % and 
an absolute calibration error between orders/modules of ~ 5 %. 

The mid-IR (MIR) SED consists of a photospheric continuum component and various 
dust features on top of it. To obtain a clear view on the signatures of the dust we subtracted 
a blackbody curve representing the photosphere. For this purpose we obtained parallel near 
infrared photometry (JHKL) at Siding Spring Observatory (SSO). The blackbody tempera- 
tures (3400-3900 K) were derived from J — K values using the relation by Houdashelt et al. 
(2000). The flux levels of the blackbodies were chosen such as to fit the near infrared flux 
levels. For sample stars with low mass loss rates we assume that the infrared excess will not 
significantly affect J — K. Fig. 2 gives an example of a spectrum before and after blackbody 
subtraction. In the following we will present and discuss only continuum subtracted spectra. 
The chosen continuum subtraction method naturally introduces some uncertainty in the de- 
rived residual dust emission, but this uncertainty does not affect the positions and relative 
strengths of the various dust features on which this paper is focused. 
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3. The MIR spectra of AGB stars in 47 Tuc 

3.1. Phenomenology of the dust spectra 

V5, V6 and V7 can be fitted nicely with single blackbodies of about 3700 to 3900 K and 
must be nearly dust-free. Dust features were detected in the remaining 8 of the 11 sample 
stars. The three dominant dust features are located at approximately 9.7, 11.5 and 13 /im. 
Additionally, broad features at 18-20 //m are visible. However, not all features are visible in 
all stars. In Fig. 1, we show the dust spectra and the location of the corresponding stars 
in the logP-K diagram. V13 is presented separately in Fig. 2. From the logP-K values we 
can form three groups among the sample stars. The spectra of pairs V5/V6, V4/V8 and 
V1/V2 are very similar, suggesting a typical spectrum is associated with a specific location 
in the logP-i^-diagram. An interesting case is V3 with similar logP-i^ values as V1/V2. 
The V3 spectrum has a very weak dust feature sitting on top of a broad mid-IR emission 
with a wavelength dependent slope. This slope can be roughly approximated by a blackbody 
curve on top of which a weak 9.7//m feature stands out. While not obvious, this is the same 
spectral signature as in V1/V2. In 13.31 we further discuss these differences as a function of 
stellar variability. 

The change of the dust features along the AGB is in agreement with earlier studies on 
the relation between dust properties and stellar parameters. While previous works on this 
subject (e.g. Onaka et al. 1989, Sloan & Price 1996, Hron et al. 1997, Heras & Hony 2005) 
refer to an inhomogeneous samples of stars, our GC spectra minimize possible effects of 
such inhomogeneities. At the lower luminosity end of the AGB, hardly any dust is present. 
Effective dust production seems to start only after the star reaches the luminosity of V21 
corresponding to a value of approximately 20001© (cf. Lebzelter & Wood 2005). The 13 fim 
feature clearly dominates the spectrum at this luminosity. Changing from V21 to V4 brings 
about a switch in pulsation mode from overtone to fundamental mode and a strong increase 
in the observed velocity amplitude. The related change in the MIR spectrum is an increase in 
the 11.5 fim feature, a decrease in the strength of the 13 /um, and the occurrence of a 9.7 /im 
feature. Going to even higher luminosities leads to VI with its dominating 9.7 /im feature 
and only a very weak indication for the 11.5 /im feature, sitting on the long wavelength side 
of the 9.7 /im emission. 

Summarizing these observations, we can say that the 13 /im feature appears in the least 
evolved stars, later looses importance relative to the 11.5 /im peak and a feature at 9.7/im 
and finally vanishes completely in the highest luminosity stars. The latter are dominated 
completely by the 9.7/tm emission. 
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3.2. Identification of the dust features 

Among the observed emission bands only the 9.7 /xm feature has a generally accepted 
identification, amorphous Mg-Fe-silicate (see, e.g., Dorschner et al. 1995 for optical con- 
stants). No consensus exists on the carriers for the broad 11.5 /xm band, the narrow 13 fim 
band, and the various features between 17 and 20 fim. Since the sample stars are all spectral 
type K or M the dust species must be oxides. 

First detected in IRAS spectra (Little-Marenin & Little 1988), the 13 //m band was 
identified as part of a 'three-component-feature' consisting of a 9.7, ~11, and the 13 /xm 
band itself. Our spectra of 47 Tuc V4 and V8 show the 'three-component-feature'. The 
13 yum feature is likely an Al-0 stretching vibration. The carrier has been identified as either 
corundum (a-Al203, Glaccum 1995, DePew, Speck & Dijkstra 2006) or spinel (MgAl204, 
Posch et al. 1999, Fabian et al. 2001, Heras & Hony 2005). 

The origin of the broad 11.5 /im band in spectra of 0-rich AGB stars has been extensively 
discussed (e.g. Stencel et al. 1990, Lorenz-Martins & Pompeia 2000, Egan & Sloan 2001, 
Posch et al. 2002, Heras & Hony 2005). Although it is especially difficult to identify uniquely 
such a broad solid state band, there remains little doubt that amorphous AI2O3 is the band 
carrier in this casej^ 

The band (or double-band) at 18-20 /im is most distinctly seen in VI (18yum), V4 
(18+20 yum), V8 (20 /xm), V13 (20 /xm), V18 (18+20 /im) and V21 (20 /im). The observed 
differences in the positions of the feature's maxima probably indicate a transition from a 
crystalline oxide (with a composition between MgO to FeO, see Posch et al. 2002) as the 
dominant opacity source in the 19-21 /im region to an amorphous silicate's bending vibration 
as the primary emission mechanism. Sloan et al. (2003) found a correlation between the 
strengths of the 13 and 20 /im features. 

In Fig. 2 we show the spectrum of V13 with the band identifications for the oxide dust 
features suggested above. The spectrum of V13 is further discussed in section 13. 4[ 

3.3. Variability of the mid-IR spectra 

At least for fundamental-mode AGB pulsators cyclic variability leads to significant 
changes in size and Tgfj of the star. These could effect the observed MIR spectra through 



^ Amorphous AI2O3 is sometimes called 'corundum' in the literature. However, this name should be 
reserved for the crystalline a-phase of AI2O3. 
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cyclic changes in both dust production and photospheric flux. From Spitzer we have only 
one observational epoch for each star. However, additional information provides indications 
of the impact on the mid-infrared spectrum. 

Our sample contains the three miras, VI, V2, and V3, all sharing a similar location in 
the logP-K-diagram. Their spectra show the same dust features, but with a very different 
contrast against the photospheric background. From the simultaneous K photometry, we 
estimate that VI was observed close to light maximum, while V2 and V3 were close to 
minima. We speculate that the background is related to the pulsation phase. This is seen in 
model atmospheres (Aringer 2005). Phase dependent background can possibly explain the 
difference between our spectrum of V8 and one published by van Loon et al. (2006). While 
dust features are quite obvious in the Spitzer observation, the van Loon spectrum taken from 
the ground shows no features at all. The K amplitude of V8 is too small to safely estimate 
the phase from the observed K magnitude. 



3.4. Two peculiar cases: V13 and V18 

Two stars do not seem to fit in the dust evolution scenario described in I3.lt V13 and 
V18 which both have IR excess emission much too strong in relation to their position in the 
logP-K diagram. The residual dust emission spectrum of V13 (Fig. 2) is almost devoid of a 
10 yum band. This alone is not surprising, because it also is the case for V21. However, V13 
is almost Imag fainter in K than V21. The spectra of both V13 and V21 are characterized 
by pronounced 11.5, 13 and 20 //m peaks. VIS shows a spectrum very similar to the mira 
VI, i.e. it is dominated by a 9.7/im silicate dust feature, yet is ~1.25 magnitude fainter at 
K. 

What properties of V13 enable the formation of dust? The only remarkable pulsational 
characteristic of V13 is its long secondary period. This long period is also clearly visible in 
the radial velocity (Lebzelter et al. 2005) as well as photometric variation. The length of the 
secondary period has not be determined accurately but seems to fall onto the Wood et al. 
(1999) P-L sequence "D". The interpretation of this sequence in terms of stellar pulsation 
or binarity is still not clear (e.g. Wood 2006). 

We interpret the dust spectrum of V18 in the light of the earlier suggestion that this star 
is currently undergoing a thermal pulse event. VI 8 would then have an interpulse luminosity 
and dust spectrum similar to that of an 47 Tuc mira, e.g. VI. According to evolution models 
of thermally pulsing AGB stars by Vassiliadis & Wood (1993), a significant disturbance of 
the typical pulsation period of a star occurs during the thermal pulse (cf. their Fig. 3), which 
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would explain the strange location of V18 in the logP-i^ diagram. 

4. Consequences for the evolution of dust on the AGB 

Our evolutionary sequence of the dust spectra indicates that oxygen-rich AGB stars first 
form dust rich in oxides of aluminum and magnesium (specifically, AI2O3-, MgAl204- and 
(Mg,Fe)0-). As a star climbs up the AGB, the original aluminum-magnesium oxide dust 
contributes less to the emergent dust shell spectra, while the relative strength of amorphous 
silicate bands, especially the one at 9.7 yum, increases. This scenario is in agreement with 
considerations published by Stencel et al. (1990) linking this evolutionary sequence to the 
higher electron affinity of Al (compared to Si) to oxygen. It is the high binding energy of solid 
Al oxides which enables them to condense at much higher temperatures than silicates (sec, 
e.g., Gail 2003, Fig. 4). Once the Al oxides have formed, silicates will start to precipitate on 
their surfaces which gradually leads to a transition to a silicate dominated dust spectrum. 
Hron et al. (1997) pointed out that the 13 /xm feature reaches significant strength only in a 
narrow range of effective temperatures and optical depths. They proposed an evolutionary 
sequence in the course of which the carrier of this band would be gradually incorporated 
into amorphous silicates with increasing mass-loss rate. Most recent model calculations by 
Woitke (A&A submitted) indicate that AI2O3 can exist closer to the star while silicates are 
formed at larger distances, i.e. once effective mass loss has set in. 

A similar scenario for aluminum-magnesium oxide/silicate dust evolution has been pro- 
posed by Onaka et al. (1989). They derive a relation between the asymmetry factor of the 
light curve and the importance of the various dust features with the stars exhibiting the 
strongest asymmetries showing the most pronounced silicate emission. The asymmetry of 
the light curve can be related to the strengths of the shock front occurring during the cyclic 
pulsation as discussed by Onaka and coworkers. Lebzelter et al. (2005) showed that the 
near-IR velocity amplitude increases along the AGB of 47 Tuc. Thus a relation between 
pulsational properties and the dust composition can be expected. The necessity of taking 
into account pulsation characteristics as an additional factor beside luminosity is illustrated 
by comparing the spectra of V21 and V4 (Fig. 1). The overall shape of the mid-infrared 
spectrum - but not the dust features that occur - seems additionally to depend on the pul- 
sational phase. A more complete discussion of the individual dust spectra will be given in 
a forthcoming paper, where we will give a more detailed presentation of the spectra and a 
comparison with combined atmospheric and dust models. 

This work is based on observations made with the Spitzer Space Telescope, which is operated 
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Fig. 1. — Development of dust spectra in 47 Tuc. Right and lower left: Spitzer spectra of 
the AGB stars in 47 Tuc (for V13 see Fig. 2) with blackbody curves subtracted (see text). 
Spectra are plotted as individual data points with errors bars derived from the difference 
between the two nod positions. The features at 9.7, 11.5 and 13 jiva. are indicated by vertical 
dashed lines. Top left corner: logP-K-diagram of the sample stars adapted from Lebzelter et 
al. (2005). Solid lines mark the locations of overtone (B) and fundamental mode (C) radial 
pulsators (Wood et al. 1999). V13 is plotted with both of its periods. 
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Fig. 2. — The Spitzer mid-IR spectrum of VI 3 (open symbols) and the residual dust emission 
(filled symbols). The grey line gives the blackbody used for continuum subtraction. Error 
bars are derived from the difference of the two nod positions. The indicated dust feature 
peak positions are from the optical constants published for amorphous AI2O3 (Eriksson et 
al. 1981), crystalline MgAl204 (Fabian et al. 2001), Mgo.1Feo.9O (Henning et al. 1995), and 
(hardly detectable) amorphous sihcate (Ossenkopf et al. 1992). 



